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This lecture is a tutorial introduction to coherent effects in disordered 
electronic systems. Avoiding technicalities as most as possible, I present some 
personal points of view to describe well-known signatures of phase coherence 
like weak localization correction or universal conductance fluctuations. I show 
how these physical properties of phase coherent conductors can be simply 
related to the classical return probability for a diffusive particle. The diffusion 
equation is then solved in various appropriate geometries and in the presence 
of a magnetic field. The important notion of quantum crossing is developed, 
which is at the origin of the quantum effects. The analogy with optics is 
exploited and the relation between universal conductance fiuctuations and 
speckle fluctuations in optics is explained. The last part concerns the effect 
of electron-electron interactions. Using the same simple description, I derive 
qualitatively the expressions of the Altshuler-Aronov anomaly of the density 
of states, and of the correction to the conductivity. The last part, slightly 
more technical, adresses the question of the lifetime of a quasiparticle in a 
disordered metal. 



1 Introduction : phase coherence and disorder 

Although the topic of this School mainly concerns nanoscopic systems, this 
set of lectures is devoted to an intermediate range, between the nanoscopic 
and macroscopic scales, the so-called mesoscopic regime [1]. In this regime, 
the system to be considered may be large compared to the mean free path of 
the electrons. Disorder plays then a very important role and, in the so-called 
diffusive regime, the interplay between disorder and quantum interference 
effects is crucial. This is the main subject of these lectures. Here, electronic 
interactions will be treated as a perturbation, in contrast with other topics 
discussed in this School where the electronic correlations may play the most 
important role. I will try to present some personal points of view in order to 
describe these well-known signatures of phase coherence like weak localization 
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or universal conductance fluctuations. The goal of these lectures is to avoid 
technicalities as most as possible. The last part concerns the effect of electron- 
electron interactions. 

To describe interference effects in electronics, it is useful to compare with 
simple facts known in optics. The simplest experiment with light is the two- 
slit Young experiment and phase coherent effects considered in electronics 
are nothing but some more sophisticated versions of the Young experiment. 
This two-slit experiment can be also performed with electrons in vacuum [2] 
but here we shall consider metallic wires, that are complex disordered media. 

In vacuum, an electron beam is split in two parts and the intensity is 
measured on a screen. The topological equivalent in a metal consists in a 
loop pierced by a magnetic flux and we measure the current resulting from 
the interferences between the two paths, see figure [1] In optics the way to 
probe the interference pattern on the screen is to change the optical path 
between the two trajectories, by changing the nature of the medium, that is 
its optical index. For electrons, the charge is coupled to the vector potential 
A and these interference pattern can be modified with a magnetic field. 



Fig. 1. ie/i ; schematic representation of the Aharonov-Bohm effect. A flux tube of 
flux 4> is placed behind the two slits. Center : Schematic equivalent of the Aharonov- 
Bohm experiment in a metallic ring. Right : Magnetoresistance oscillations of a Au 
ring and its Fourier transform 

If we try to transpose what is known from optics to electronics, we measure 
a current intensity which is proportional to the probability for the electrons 
to traverse the loop. To calculate this probability in quantum mechanics, we 
have to add the contributions of two quantum amplitudes corresponding to 
the two sides of the loop, and the current (a probability) is proportional to 
the square of this quantum amplitude. Each quantum amplitude has a 
phase (pi : 
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We have to sum the amplitudes and take the modulus square. For the current, 
we get a classical term plus an interference term : 

/cx|V'i+V'2|' = 2V''[l + cos(0i-02)] . (2) 

Classically, the current should be the sum of the two currents, this cor- 
responds to Ohm's law, but quantum mechanically there is some addi- 
tional phase effect. And if we apply a magnetic field (this is the well-known 
Aharonov-Bohm effect [5]), there is an additional phase along each of the two 
trajectories 

dipi^^J^A-dl , 6ip2^^J^A-dl , (3) 

so that the phase difference is modulated by the circulation of the vector 
potential along the closed loop formed by the two trajectories, 

A(p = Sipi ~ S(p2 = T (f A - dl = 27r-f- , (4) 

and is proportional to the magnetic flux through the loop. 0o = h/e is 
the fiux quantum. This tells us that the wave functions, the energy levels, 
thermodynamic and transport properties must be function of this flux with 
period (po- Of course this does not tell us anything about the amplitude of 
these oscillations. Unlike the case of the Young experiment with light or for 
electrons in vacuum, the problem here is much more complicated due to 
disorder. The question now is : do interference effect survive in the presence 
of disorder? Do oscillations persist? 

Indeed, some oscillations remain in the presence of disorder. The pionner 
experiment founding the field of mesoscopic physics was performed by Webb 
et al. in 1985 [3]. They measured the resistance of a ring of micronic size, in 
the presence of a magnetic field (figure [T]). They found that the resistance 
oscillates with the field, proving the existence of an interference effect, even in 
the presence of disorder. Two interesting features must be noted : firstly, the 
period of the oscillations is the fiux quantum 00 as expected. Secondly, the 
typical amplitude of the oscillations of the conductance is AG — AR/R^ ~ 

The necessary condition for these oscillations to exist is that phase co- 
herence is preserved, a condition which is obeyed in vacuum. But in a metal, 
this phase coherence is broken because electrons interact with other degrees 
of freedom {e.g. phonons, other electrons, magnetic impurities). Because of 
this coupling, each electron can lose the memory of its phase. This happens 
on a typical length, called the phase coherence length, denoted by L^, which 
depends on the coupling to these degrees of freedom. Typically at IK, it is 
of order of 1 /i. This is the mesoscopic range. Of course, such oscillations 
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do not exist for a macroscopic ring, because phase coherence is fost at the 
macroscopic level. 

There is another important scale, due to static disorder, over which elec- 
trons experience collisions over impurities. This scale le is named the elastic 
mean free path. It is much smaller than the phase coherence length. Since elas- 
tic collisions do not break phase coherence, electrons may experience many 
collisions without losing the memory of their phase. Each electron keeps the 
memory of its phase typically on a scale le- The physics we are going 

to discuss corresponds to length scales which are much larger than the elastic 
mean free path l^., but smaller than L^, so that the sample can be consid- 
ered as quantum mechanically coherent. This regime is called the mesoscopic 
regime. 

We want to construct a theory for these oscillations. How to explain their 
amplitude? How can they survive disorder? Another important key exper- 
iment was performed by Sharvin and Sharvin on a cylinder [4], before the 
Webb's experiment on a single ring. In this case, there are also oscillations, 
but instead of being (/fg-periodic, they have a period ^o/2. The fundamen- 
tal (f>o disappears and the experiment reveals the second harmonics of these 
oscillations. 

There is a simple way to understand this frequency doubling : assume 
that the cylinder can be viewed as a pile of independent rings, it realizes an 
average of the oscillations of several independent rings (similar effect for a 
network of rings [6]). For a given ring, the oscillations have a given phase 
which depends on the disorder configuration (This phase is or tt for a two- 
terminal geometry [7]). For another ring, the oscillations have a different 
phase. When averaging over several rings (as is done in a cylinder), because 
of this random phase, the oscillations vanish. 

So, disorder is destructive for quantum interferences. However this simple 
argument would tell us that no oscillations should survive disorder averaging. 
This is not the case since the cylinder experiment shows oscillations with 
period (/'o/2. In average, there are still oscillations, but with period h/2e. 
This simple and very important fact tells us that some contributions survive 
disorder averaging. So the question is : how is it possible that some robust 
contribution survive disorder averaging ? 



2 Important scales 

We shall consider weakly disordered metals, such that the average distance 
between two collision events is much larger than the Fermi wave length : 1^ ^ 
Xp. This condition allows for a semiclassical description of electronic waves. 
Moreover we assume that the typical size L of the system is much larger than 
the mean free path l,, so that the electronic motion in the sample is diffusive : 
electrons collide elastically many times while traversing the system. Finally 
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Webb et al. experiment 


Sharvin-Sharvin experiment 






Even in the presence of disor- 
der, phase coherence is preserved 
on distances much larger than 
— > sample specific interference 
effects, period 4>o- 


Some interference effects survive 
disorder averaging, period 0o/2 
— > pairing of trajectories. 



Fig. 2. Conclusions that can be drawn from the two pionniering experiments in 
mesoscopic physics. 

we assume that the system is completely phase coherent, that is L ^ L^. To 
summarize, the system is weakly disordered, diffusive and mesoscopic : 

Ai^ < ?e < L < . (5) 

We shall see that phase coherent effects may not disappear but are simply 
reduced in the macroscopic limit < L, so that we shall also consider the 
case : 

Ai^ < ?e < ^0 < i ■ (6) 

In strong disorder, when the mean free path becomes of order of the Fermi 
wave length kplf, ~ 1, interference effects are strong and lead to localization 
of the electronic waves. This is the domain where the electronic states are 
exponentially localized in space, with Anderson localization from extended to 
localized waves. This topic will not be covered here. In the opposite regime of 
very weak disorder, the mean free path becomes so large (or the system is so 
small) that the mean free path becomes larger that the system size. Collisions 
with impurities are rare, and occur mainly on the boundaries of the system. 
This the so-called ballistic regime, where the physics is mainly driven by the 
structure of the boundaries, i.e. the shape of the system. For most shapes, the 
trajectories are chaotic. A common method to describe this regime is the so- 
called Random Matrix Theory also used in other fields of physics like nuclear 
physics. There are quite interesting common features between some aspects 
of transport in chaotic dots and nuclear physics. Let us also emphasize that 
the Random Matrix Theory of scattering or transmission matrices can also 
be used to describe diffusive systems [8]. 

Since the typical size L of the system is much larger than the mean free 
path, the electronic motion is diffusive. The average distance between collision 
events, the mean free path is related to the collision time Te, le = vpTe, vp 
being the Fermi velocity, since the electronic motion between two collisions 
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is ballistic. For times much longer than the collision time Tg, the motion is 
diffusive and the typical distance an electron can reach after a time t scales 
like 

r^^Dt , (7) 

where D is the diffusion coefficient given by D — vple/d, d being the space 
dimensionality. This relation tells us that for a finite system of size L, a very 
important scale appears : this is the time for which an electron typically sees 
the boundaries of the system. It is called the traversal time, or Thouless time. 
It is the time for an electron to "realize" that the system is finite. It is given 

by 

rn = ^ . (8) 

To this characteristic time, is associated a characteristic energy, the Thouless 
energy Ec : 

h hD 

Ec = — ^-rj. 9 
td 

This energy scale plays a major role in the description of thermodynamic and 
transport properties of mesoscopic diffusive systems. For time scales smaller 
than Td, the electron propagates like in infinite space. The diffusive motion 
depends on the space dimensionality of the system. On the other hand, in 
the long time scale, the electronic motion explores the entire system, this is 
the so-called ergodic regime. 

Using eq.(l7]), we can associate to the phase coherence length a char- 
acteristic time, the phase coherence time : 

r, = § . (10) 
This is the time during which an electron keeps the memory of its phase. 



K agccdc >K ftee difflisbn X — baIList±: 



A fi/i; E^ 

Fig. 3. Characteristic energy scales defining the different regimes studied in coher- 
ent multiple scattering. Explanations are given in the text. 



Figure [3] presents a scale of characteristic energies (or inverse characteris- 
tic times). At short time scales, the motion is ballistic. For times larger that 
Te, the motion is diffusive in free space. Above r^i, the motion is bounded, 
this is the ergodic regime. Then separates the mesoscopic regime and the 
classical regime. If > td {L^ > L), the system is mesoscopic. 
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In the diffusive regime, the space dimensionality d plays an important role. 
Moreover, the one-dimensional case is somehow special. In strictly one dimen- 
sion, it is known that there no diffusive regime since, in the presence of disor- 
der, all states are exponentially localized. So we shall not consider this case 
but rather gwasi-one-dimensional systems, with a transverse width a, so that 
the real motion is three-dimensional (the proper Schrodinger equation to be 
solved would be three-dimensional with a quantization of the transverse com- 
ponent of the wave vector), but the diffusion is one-dimensional. Instead of 
having one transverse propagation channel, there are many transverse chan- 
nels (the transverse size is much larger than Xp). At short time (smaller than 
a "transverse" Thouless time a^/D), diffusion is three-dimensional, but at 
larger times, there is a one-dimensional propagation of the diffusion cloud. 
When considering the transport through a wire (L 3> a), we shall be inter- 
ested in time scales necessary to traverse the wire, that is times of order of 
td, which is much larger than the transverse time, so that we can consider 
that at this time scale the diffusion is one-dimensional. 



3 Classical probability and diffusion equation 

The aim of these lectures is to propose a qualitative description of physical 
phenomena, trying to avoid sophisticated tools and keeping in mind that we 
are essentially concerned by the calculation of disordered averaged quantities. 

An average quantity like the conductance basically measures the proba- 
bility for electrons to cross the system. What is the nature of this probability? 
Let us first spend some time to describe the probability P{r,r') which de- 
scribes the propagation of a particle from a point r to a point r' . In quantum 
mechanics, this propagation is described by a probability amplitude. This 
amplitude is called a Green's function G{r,r'). We do not aim to develop 
the theory of Green's functions. For our purpose here, it is sufficient to note 
that there are many possible scattering trajectories from r to r' . Thus a 
Green's function has the following structure : it is the sum of all the prob- 
ability amplitudes corresponding to various multiple scattering trajectories 
from r to r', each trajectory being characterized by an amplitude and a phase 
proportional to its action, that is its length |9] : 



Now, we want to know the probability to find a particle at point r' if it has 
been injected at point r. The probability to go from r to r' is the modulus 
square of the amplitude. From eq. pT|) . we see that this probability is the 
sum of amplitude squared terms, plus interference terms which pair different 
trajectories j and f : 




(11) 
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|G(r,rOP=^^,(r,r')A,-Kr,r')=5]|A,(r,r')P + 5]A,(r,r')A*,(r,r') 

(12) 

which an obvious generahzation of eq. ^ for the two-sht configuration. Since 
we know that, in quantum mechanics, one must add ampHtudes instead of 
intensities, the interference term (the second term in eq. I12p cannot be a 
priori neglected. This second term describes interferences between different 
trajectories j and j'. Each contribution in this sum has a random phase 
which depends on the detail of the impurity configuration. Since the phases 
are uncorrelated, at the first level of approximation, we may expect that 
the contribution of the interference term cancels upon disorder averaging. So 
quantum effects seem not to be so important because of the vanishing of this 
contribution. We shall see however that this is not exactly the case. Within 
this approximation, the second term cancels and the probability is essentially 
given by the sum of intensities : 

|G(r,r')P = ^|A,(r,r')P • (13) 

3 

We see that the phases have disappeared. So the remaining term is completely 
classical. Indeed, let us assume that some event changes the phase of the 
amplitude Aj. The complex amplitude A* gets the opposite phase, leaving 
the probability unchanged. 




r 



&(^A'} PotK'"') 

Fig. 4. Schematic representations of a Green's function G{r ,r') and of the classical 
probability Pct{r,r') oc J^^, \Aj{r,r')\^. The upper diagrams exhibit a few collision 
events, which are not represented on the lower diagrams. 

To have a simple picture of eq. [13| , we represent on figure 2] a quantum 
amplitude as a line (it is rather a sort of Brownian trajectory). Its complex 
conjugate is represented as a dashed line. The first term in eq. (jl3[) corre- 
sponds to the pairing of a trajectory with its complex conjugate, and we see 
immediately why the phase disappears. The quantity J2j ''"'W ressem- 

bles the classical probability. We call it a " Diffuson" . To be more precise, but 
without any proof, we define the probability P{r^r' ,uj) as 



P(r,r',c.) = -J_G,(r,r')G*„„(r',r) 



(14) 
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The Green's function and its complex conjugate are taken at different energies 
(or frequencies) e and e — w. One can check that this probabihty is correctly 
normalized, that is P{r,r' ,t)dr' = 1, where P{r,r',t) is the Fourier 
transform of P{r, r', lo). Starting from the Schrodinger equation in a random 
potential and after disorder averaging, it is possible to show that in the limit 
kple ^ 1, the probability P{r,r',iv) defined by (fT4|) reduces to the Diffuson 
Pci{r,r') cx |Aj(r, r')p. For slow spatial variations, Pci{r,r' ,uj) is the 
solution of a classical diffusion equation : 

{-iuj~DA)P,i{r,r',uj) = 5ir-r') (15) 

where D is the diffusion coefhcient. Doing this, we have only considered 
classical contributions to the average (HH). We shall study later the corrections 
to this classical probability. 

Among the solutions of this diffusion equation, one very important is 
the return probability which enters in many physical quantities. It is the 
probability P{r^r,t) for an electron to return to its original position after 
time t. I will also consider the space integrated return probability : 

P{t) = J P{r,r,t)dr . (16) 



In free space, the solutions of equation (1151) simply obtained from the 
Fourier transform 

l + Dq^'^P{q,t)=Sit) (17) 
whose solution P{q,t) is simply 

P(q,t) = e-^*'' . (18) 
Fourier transforming back, we find easily 

so that the return probability is given by 

P(r,r,i) = - 5—— and P(t) = -, ^^-77^ . (20) 

where f2 is the volume of the system. The dependence on the dimensionality 
d of the return probability is crucial since it will explain why dimensionality 
plays a so important role in mesoscopic physics of diffusive systems. 



4 Conductance 



Now I wish to come to very simple and qualitative considerations about the 
conductance of a disordered system, which will be useful for these lectures. 
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4.1 Classical conductance as the ratio of two volumes 

Consider the conductance G. Since it has the dimensions of e^//i, we can 
introduce a dimensionless conductance g as 

g = GI{se'lh) . (21) 

Since this quantity is dimensionless, it may be written as the ratio of two 
physical quantities. For example, by simple manipulations, it can be written 
as the ratio of two energies : g cx E^j the Thouless energy and the 
average level spacing A. Here I would like to write it as the ratio of two 
volumes. Let us start with the classical Drude conductivity uo- From Einstein 
relation, it is given by 

ao = se^DpQ , (22) 

where D is the diffusion coefficient and po is the density of states at the Fermi 
level for one spin direction. The factor s = 2 accounts for spin degeneracy. 
By Ohm's law, the conductance G for a three-dimensional system is given 
by G = aoS/L, S being the section and L the length of the sample. More 
generally for a hypercube of typical size L in d dimensions, it is given by 
G = aoL'^~^. Introducing the Thouless time defined by (O, let us rewrite 
the conductance as 



G^se''poL''/TD . (23) 

The density of states at the Fermi level po can be written as po = dAd/2'KX'^^hvF, 
where Aj^ is the Fermi wavelength, vp is the Fermi velocity, and is the 
volume of the unit sphere {A3 = 47r/3, A2 — tt, Ai — 2). An easy way to 
recover immediately this result is to say that the total number of states is 
{kpLY, so that by derivation with respect to the energy, we have necessar- 
ily Po cx kp/ep — klpT^ /hvp oc l/hypX^^ . As a result, the dimensionless 
conductance g can be written as 

9 = d^d , (24) 

where J7 = i'' is the volume of the system. The dimensionless conductance 
quantity appears as the ratio of two volumes, the volume f2 of the system 
and the volume of a tube of length vpTn and of section A^"""^. We shall see 
later that this formulation will be quite useful to measure the importance of 
interference effects. 



4.2 Conductance and transmission 

Our starting point to describe electric transport is the Landauer formalism. 
Even staying at a very qualitative level, this formalism is quite natural since 
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Fig. 5. In the Landauer formalism, the conductance is related to the transmission 
coefficient between different incoming and outgoing channels. 



it expresses the conductance as a transmission coefficient through the disor- 
dered sample. 

Consider a disordered conductor of length L and section S = W''"'^ . 
It is connected to perfect conductors (figure [5]) which can be considered as 
wave guides where free electronic waves propagate. In this geometry, the 
transverse wave vectors of the eigenmodes (also called channels) are quantized 
by transverse boundary conditions. One can define a transmission coefficient 
Tab from an incoming channel a (ingoing wave vector fc^) to an outgoing 
channel b (wave vector fe^). The Landauer formula reads : 

G = s^^Ta6 . (25) 

a,b 

To calculate the number of transverse channels, one considers that electrons 
are injected at the Fermi energy, i.e. such that \ka\ = |fch| = kp. The trans- 
verse component is quantized in units of 2tt/W. This quantization imposes 
the number of channels. In d = 2 and d = 3, their number is 

27rfcf? Trfcp k1,S , , 



Let us consider now the structure of the transmission coefficient T^h- It 
is the square of an amplitude and it has, with minor differences, the same 
structure as the probability P{r,r',uj). The main difference is the following : 
instead of injecting a particle at a point r inside the sample, we inject a 
plane wave ka from outside the sample. In particular, the boundary condi- 
tions have to be treated properly. But, without entering into details, we may 
easily understand that, after disorder averaging, the average transmission co- 
efficient and consequently the conductance can be related to the probability 
to cross the sample. More precisely for a 3d sample, one can show that the 
dimensionless conductance is [T] 



g = -MvFP{0,L) 



(27) 
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where P(0, L) is the solution of the diffusion equation (fT5|) with appropriate 
boundary conditions. It is given by P{0, L) — l^/DL so that 



To obtain these results quantitatively, there are some technicalities that we do 
not describe here [T] . What should be remembered is the message of fig-IH]: the 
conductance is proportional to the classical probability to cross the sample. 
This statement is sufficient to understand how coherence effects appear. 



Fig. 6. The conductance is proportional to the classical probability to transmit 
channel a to channel b (summed over channels). The object which represents this 
probability is the sum of contributions of paired trajectories as introduced in figure 
g) We call It a "Diffuson". 



5 Quantum corrections and quantum crossings 

The Diffuson is a classical object. It does not depend on the phases of the 
complex amplitudes. In the diffusive regime, it is solution of a diffusion equa- 
tion. However, we may have to check whether we have not left aside additional 
effects when throwing out all the interference terms in relation p^ . It turns 
out that some of these terms have quite interesting consequences. 



Fig. 7. a) Crossing of two Diffusons. b) Detail : the volume of the intersection 



Indeed, quantum effets can appear when two Diffusons cross, or when 
a Diffuson crosses with itself. The notion of quantum crossing is extremely 
important because it is the source of quantum effects. The Diffuson being a 




(28) 





region is proportional to ^le- 
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classical object, coherence effects can only appear because of these quantum 
crossings. They are at the origin of the weak localization correction and of 
universal conductance fluctuations. Let us try to get some intuition about 
these crossing events. 

Figure [7la shows that a crossing mixes four complex amplitudes which 
belong to two incoming Diffusons and pair them differently. The two emerging 
Diffusons are built with amplitudes Aj and Aji coming respectively from each 
of the incoming Diffusons. They have the same phase since they follow the 
same path. The quantum crossing, also often called a Hikami box in a more 
technical context, is thus an object whose role is to permute the quantum 
amplitudes. It is necessarily short ranged, because trajectories have to be 
as close as possible to each other to avoid dephasing (figure [7lb). Since it 
appears between two successive collisions on impurities, and since the phase 
mismatch between trajectories has to be smaller that 27r, one sees that the 
volume of this object is of order X'^^le- 

It is important to evaluate the probability of occurrence of such quantum 
crossings because it will be shown to be a measure of the importance of 
quantum effects. Since the volume of a quantum crossing is of order X'^^lg, 
a Diffuson propagating during a time t can be seen as an effective object 
of length C = vpt and of "cross section" A^~^. Thus, it has a finite volume 
vpy^^t. The probability dpx{t) of crossing of two Diffusons after a time dt 
in a volume 12 = L'' is thus proportional to the ratio between the volume of 
a Diffuson and the volume of the system : 

Xir^vpdt 1 dt 

dpx (t) = ^ „ (X 29 

S2 g TD 

where we have used (|24l) to introduce the dimensionless conductance g. Con- 
sider now an open system coupled to reservoirs. The time needed to travel 
throughout the sample is the Thouless time tu — L^/D. The probability of 
quantum crossing during this time is given by 

Px{td)^ dpx{t) ^ „ =i - . (30) 



Jo ^9 
This is exactly the inverse conductance ! I believe that this is the most im- 
portant message to understand phase coherence effects in disordered systems. 
All these effects can be simply understood in terms of quantum crossings and 
the probability of such crossings which measures the importance of quantum 
mechanical effects is simply given by the inverse of the dimensionless conduc- 
tance g. 

In a good metal, the conductance g is large, the volume of the tube is 
small, electrons do not spend much time in the system and quantum effects 
are very small. In the opposite limit, when g becomes of order 1, the volume 
of this tube is of the order of the volume of the system. It is so big that 
the probability of quantum crossing is of order 1. This corresponds to the 
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Anderson regime where electronic waves are localized by strong disorder. 
Here we shall not consider this regime but only the small disorder regime 
where quantum effects remain small. The approach to Anderson localization 
can be viewed as the proliferation of quantum crossings. 

As a first qualitative but important conclusion of our discussion, we see 
that classical transport is described by a conductance G — sge^/h. Quan- 
tum corrections are smaller than classical terms by a ratio 1/g. This tells us 
immediately that the quantum corrections are of order G/g, that is e^/h ! 

6 Weak localization 

6.1 Weak localization and quantum crossings 

We have seen that the classical probability and the conductance can be ex- 
pressed as a sum of contributions of pairs of complex conjugated trajectories. 
Since trajectories can have quantum crossings, they can form closed loops 
(figure [8]). It turns out that in such a loop (whose contribution is not in- 
cluded in eg. [T5)) . the trajectories are time-reversed. One trajectory j and its 
time reversed go in opposite directions. If there is time-reversal symmetry, 
they have the same action and thus they have exactly the same phase. This 
phase can be quite complicated because it depends on the disorder configu- 
ration but it is the same for both trajectories. So the contribution of these 
loops does not cancel on average. If the end points are far away like in fig. [51 
the contribution of these new trajectories is small, of order 1/g, but it leads 
to an experimentally observable eff'ect : the weak localization correction to 
the conductance. This is a phase coherent effect because only trajectories of 
size smaller than the phase coherence length contribute to this additional 
contribution. 

At this point, I want to stress that many presentations of weak localization 
correction emphasize the existence of a loop of opposite trajectories, but do 
not insist on the structure of the quantum crossing. This is rather quite 
important, because this is what explains the amplitude 1/g of the correction. 
This is where phase coherence is lost. 



ft. 

ft- 




Fig. 8. Trajectory with, a quantum crossing and a loop. In ttie loop, tfie two propa- 
gations are time reversed. 



Coherence and interactions in diffusive systems 



15 



Using the same type of argument as in the previous section, let us evalu- 
ate the probability to have a loop for a trajectory which travels through the 
sample. Since there is a quantum crossing, the probability is small, of order 
1/g. Moreover, it depends on the distribution of loops in the disordered sys- 
tem. Let us call it Pint{t)- For the probability of traversing the sample with 
a loop, we have : 

1 dt 

Po{td)^ Pzntit) dp,Xt) ^ - i^„,(t)_ . (31) 
Jo 9 Jo TD 

We have also to remember that, because of decoherence in the loops, only 
those with time t smaller than contribute. The resulting probability to 
have trajectories with loops of time smaller than is : 



Po{t^)= P^nt{t) dp^{t) = - Pint{t)— (32) 

Jo 9 Jo TD 

where Pint [t) is the probability to have loops of time t. This leads to a relative 
correction to the conductivity (or to the conductance) given by 

= 7=;- = -Po{t4>) (33) 

The sign of the correction is negative because the trajectories j and j'^ have 
opposite momenta. This quantum correction to the classical Drude conduc- 
tivity is called the weak localization correction [10l|Tll ?]. The phase coherence 
is broken by the coupling of the electrons to other degrees of freedom or due 
to electron-electron interactions. Such coherence breakdown is temperature 
dependent and can be phenomenologically described by a temperature de- 
pendent phase coherence length ^^(T) = ^J Dt^{T) : trajectories larger than 
do not contribute to the weak localization correction. 

As we have seen above, the amplitude of the correction is proportional 
to Pint{t), the distribution of loops. This quantity is nothing but the return 
probability after a time t. It is not exactly the classical return probability, 
which is the product of an amplitude with its complex conjugate (Diffu- 
son). Here it corresponds to the product of an amplitude with the complex 
conjugate time-reversed amplitude. This object is named a " Cooperon" . For 
closed trajectories, and with time-reversal symmetry, the two contributions, 
Cooperon and Diffuson, are equal. The return probability is thus doubled 
due to quantum coherence. 

Eqs. ([32j[33|) have a meaning only in the diffusive regime for which t > Te 
(otherwise a loop cannot be formed). The contribution of the return probabil- 
ity has thus to be integrated between Te, the smallest time for diffusion, and 
the phase coherence time r^. Replacing the bounds by exponential cutoffs, 
the weak localization correction can be cast in the form [13] : 
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AG = -2s- P„,(t) (e-*/-*-e-*/^=)— . (34) 
n- Jo TD 

In order to evaluate AG in various situations, we now study the diffusion 
equation and its solutions. 



6.2 How to solve the diffusion equation 

In order to perform practical calculations, we have to calculate this distribu- 
tion of closed trajectories, that is, the return probability. We have to solve a 
diffusion equation for this interference part. It looks very much like a classical 
diffusion equation (I15|) . but there is an important difference. To account for 
magnetic field effects, it should be replaced by (in time representation) : 



P{r,r\t)^S{r~r')S{t) . (35) 



The effect of the field is described by a covariant derivative (like in a 
Schrodinger equation) with an effective charge 2e to account for the pair- 
ing of trajectories. 

To solve this equation (|^ . let us notice that it is a Green's equation, 
whose solutions are 

Pir,r',t)^e{t)J2rn{r)Mr')e-''-' , (36) 

n 

where 9{t) is the step function and {En,tj}n} are the eigenvalues and eigen- 
functions of the eigenvalue equation associated to ((35|) : 



9-1 p A 

-D{Wr + ^)'Mr)=EMr) ■ (37) 
nc 

From eq. (|36p. we find that the integrated probability P{t) has the simple 
form : 

P(i) = 0(t)^e-^"* . (38) 

n 

This important result tells us that in order to evaluate the weak localization 
correction in any geometry, we simply need the eigenvalues of the diffusion 
equation in the corresponding geometry. We consider now a few examples. 



6.3 Dimension dependence of the weak localization correction 

Consider an infinite system, or with size L 3> i^. For free diffusion in infinite 
space, the eigenvalues £"„ of the diffusion equation are Dq^ and the return 
probability P{t) is given by ((20| . Since P{t) is dimension dependent, we 
see than this weak localization correction depends dramatically on the space 
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dimensionality. Inserting (|20[) in 
well-known results : 



La 



, and writing = L^/D, we obtain the 



-(- 



quasi — Id 
d = 2 



(39) 



L 



Since L^(T) varies as a power-law with temperature, we obtain in particular 
the famous In T dependence of the weak localization correction in 2d. 

It show be noticed that these results are meaningful if the correction stays 
smaller than the classical conductance (which can be written in the general 
form g = Ad{kFL/2TTY^^ le/L). This defines a characteristic length ^ given 
by ^5(C) — 5; for which the weak-disorder perturbative regime breaks down. 
In Id and 2(i, it is given by 



Id 



2L 



2d 



(40) 



and for a quasi-ld system S.qid — Mlf, where M is the number of channels. ^ 
is the localization length. For a review on the strong localization regime, see 
for example [Tl] . 



6.4 Finite systems, boundary conditions 

In a mesoscopic system, the cutoff time in (l32t is provided by tu- In other 
words, the cutoff length in eqs. ([55]) is now the size L of the system instead 
of L^. Therefore, from ([55]) . we see that in quasi-ld, the weak localization 
correction is universal in the sense that it is a number, independent of disorder 
strength {le)- In 2d and Sd, the integral ([M]) diverges at small time and is cut 
off by Te, so that the correction cannot be universal [l]. 

In order to calculate quantitatively the weak localization correction and 
the return probability in a finite system, one must be careful to account 
properly for correct boundary conditions. 

If the system is closed, electrons stay inside the system, so that P{t) — > 1. 

t — *oo 

The correct boundary condition is that the probability current vanishes at the 
boundary (Neumann condition). Therefore q = nn / L with n = 0, 1, 2, 3, • • • . 
For t — > oo, the contribution of the zero mode in eq. (I38|) gives correctly 
P{t) ^ 1. 

If the sample is perfectly connected to leads, electrons can leave the sample 
and P{t) — > 0. The probability at the boundary has to vanish (Dirichlet 

t — >oo 

boundary condition) , because if it goes in the leads it never comes back in the 
same state. The zero mode in now excluded, q = rm jL with n = 1, 2, 3, • • • . 
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P{t) when t — > since the particle leaves the box at large time. Inserting 
the expression ([55]) of the return probability with i?„ = ri^Ec into gives 
immediately, in the limit — s- cx) : 



This result is proper to the perfectly connected wire. 
6.5 Magnetic field effects 

Ring or cylinder geometry : Sharvin-Sharvin oscillations 

Consider first the geometry of a ring pierced by a Aharonov-Bohm flux. In 
the presence of the flux, each closed trajectory accumulates an Aharonov- 
Bohm phase 2T:(j)/(f)Q, where is the flux through the ring. The time-reversed 
trajectory accumulates an opposite phase — 27r(/)/(/)o, so that the relative phase 
shift between the two trajectories is 47r0/0o, The fact that this relative phase 
between the two time-reversed trajectories is twice the phase enclosed by one 
trajectory is the reason why average quantities oscillate with period (/)o/2 = 



We need to calculate the return probability in this geometry. This can 
be done directly by solving ([35]). Here let us proceed by simple arguments. 
Remember that in a Id infinite space, the probability to go from r to r' 
is given by (|19p with d = I. The return probability is obtained by writing 
r — r' . On a ring, this would be the return probability without making a loop, 
1 / \/ AirDt. The return probability after one loop of perimeter L necessarily 
contains a term /iot^ rpj^g accumulated phase is 47r0/0o so that the 
flux dependent contribution of trajectories making one loop is modulated by 
cos47r(/)/0o- Adding together the contributions of m loops, we get simply the 
Fourier expansion of the flux dependent return probability : 



Each harmonics of this expansion represents the return probability after m 
loops around the ring. Inserting this expression into (|34|) and after time in- 
tegration, we obtain easily 




(41) 



h/2e. 




(42) 




(43) 



which can be resumed to obtain 



ziG(0) 



= —s 



sinh L/L^ 



(44) 



h L cosh L / — cos 4:TT(p/<j>o 



Coherence and interactions in diffusive systems 



19 



The harmonics decay exponentially with their order, since they correspond 
to longer and longer diffusive trajectories. 

For a cylinder, there is a possibility for the electrons to diffuse along the z 
axis of the cylinder, so that is simply multiplied by L^/ y/AnDt. Inserting 
this new probability into eq. ([34]), we obtain 



AGicp) = -s 



Trh L, 



In^ + 2 ^ Ko{mL/L^) cos At: m4> / 4>o 



+ 00 



(45) 



where Kq is a modified Bessel function [15]. The m = term is the usual 
2d result pQ]) . The contributions of the harmonics decay as e"^/^*. These 
oscillations, predicted by Altshuler, Aronov and Spivak, where observed by 
Sharvin and Sharvin in 1981 [illTB]. 



2d gas in a magnetic field 

In the ring geometry, all pairs of diffusive trajectories would pick the same 
phase 47r0/(^o- In a uniform magnetic field, small and large trajectories accu- 
mulate different fluxes 0(-4) = BA depending on their area A. So the return 
probability is balanced by the average {cos A'K(j}{A) / (jio) a on the distribution 
of areas A formed by the time-reversed diffusive trajectories. 

Let us start with a qualitative evaluation. Short trajectories accumulate 
a flux smaller than the flux quantum and their contribution survives. Large 
trajectories accumulate flux larger than 0o and their contribution vanishes. 
When the magnetic fleld increases, the contribution of smaller and smaller 
trajectories is progressively suppressed. Trajectories smaller than some field 
dependent characteristic length Lb corresponding to BL^g ~ ipQ will not 
contribute. To this length L b corresponds a characteristic time tb = L\ /£> ~ 
cfio/BD, so that we can expect 

(cos47r0/<?!)o)^ ~ e-*/^« . (46) 

Trajectories which enclose more than one fltix quantum do not contribute to 
the return probability. Because of this new cutoff time, we can expect a field 
dependence of the weak localization of the form 

^^^_l^^^min(L^ (47) 

TT le 

instead of 

The exact calculation is straightforward starting from eq. ([55)1 . The eigen- 
values En are solutions of an effective Schrodinger for a free particle of mass 
m = h/2D and charge — 2e in a uniform fleld B. They are precisely the 
Landau levels, namely 
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^n-("+2)— ^ ' 

where n is an integer. The degeneracy of these levels for an area 5 is 5^ = 
2^5*. The integrated return probability Pint{t) is just given by the sum 

EnSne"-^"*, that is : 

PUt,B)= . (49) 
smh(4:Tr BDt/cpo) 

where (j)Q — h/e is the flux quantum. This expression is nothing but the par- 
tition function of the harmonic oscillator. In the limit i? — > 0, one recovers 
the result for free diffusion : 5/(47rZ)<). For large times, Pint{t,B) decreases 
exponentially with the characteristic time tb — (j)o/4:T:BD introduced qual- 
itatively in ((46l) . It describes the dephasing of time reversed trajectories. 
Inserting eq. (j49p in p4p . and performing the integral, we get : 



MB) . -1; 



,1 h \ (I 



(50) 



2 AeDBTeJ \2 AeDBr^^ 

where ^ is the digamma function. This expression corresponds to the approx- 
imation (j47p . The weak localization correction is negative and cancelled by 
the magnetic field. As a result, a negative magnetoresistance is a well-known 
signature of weak localization (figure [9]). A magnetoresistance measurement 
is a very interesting and useful tool to estimate r^. The correction cancels 
when Tb — T^, that is for a field B^ corresponding to a flux quantum through 
an area B^L^. Doing the same measurement at different temperatures is the 
usual method to extract t^{T). 



7 Conductance fluctuations 

7.1 Universality as a signature of quantum coherence 

At a scale L < L^, a conductor is a quantum object. Its conductance depends 
on the interference pattern between all diffusive trajectories. This interfer- 
ence pattern can be modulated by external parameters, like a magnetic field 
or a gate voltage. For example, figure [TOla represents the variation of the 
conductance with the magnetic field, performed for 46 different samples (ac- 
tually the same sample which has been annealed several times, so that the 
impurity configuration has changed and the interference pattern is different). 
It exhibits "fluctuations" which are reproducible for a given configuration of 
disorder. They are a fingerprinf of this configuration. Figure fTOlb displays 
the average conductance, obtained by an average over the 46 samples. One 
clearly sees the weak localization correction, which is destroyed around some 
characteristic field B^. Interestingly, for the same characteristic field, the 
variance of the fluctuations, displayed on figure ITOlc is reduced by a factor 2. 
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Fig. 9. Magnetic field dependence of the magnetoresistance of a Mg film, for dif- 
ferent temperatures. The points are experimental results and the solid curves corre- 
spond to it50l) . The time t^{T) is a fitting parameter fITf . 



Universality of the conductance fluctuations is a signature of quantum 
transport. Classically, one would expect the system to be considered as an 
addition of large number of incoherent elements. This number is of order J\f = 
{L/LcY, where would be a correlation length, of the order of the mean 
free path. Then one would expect relative fluctuations of the conductance G 
of order of 



5G_ 
G 



1 



7^ 



so that, since G varies as L' 



d-2 



d/2 



(51) 



6G cx L— 



(52) 



and vanishes for large L. The system is said to be self-averaging. But the 
fact that the fluctuation 5G stays actually finite means that there are strong 
correlations due to quantum coherence. Moreover, if one considers fully co- 
herent {L < L^) systems with quite different conductances, a good metal, a 
bad metal, or a semiconductor, one finds that the amplitude of the "oscilla- 
tions" is always the same : it does not depend on the disorder. It is universal. 
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Fig. 10. Reproducible fluctuations of the 
magnetoconductance in units ofe^/h, at 
T = 45mK for Si doped GaAs. Fig. a 
shows 46 plots as function of the mag- 
netic field, for the same sample after suc- 
cessive annealing. Each plot corresponds 
to a disorder configuration and is called a 
magnetofingerprint. The amplitude of the 
fluctuations is smaller than e'^/h because 
L > (see eg. \60il . Fig. (b) presents 
the average conductance versus field. The 
weak localization correction disappears be- 
yond a characteristic field. Above the 
same field, the variance of the fluctuations 
displayed on fig. (c) is divided by a factor 
2, corresponding to the destruction of the 
Cooperon \17^ . 



BflCTeslas) 



of order /h. A priori., we are not so surprised that these fluctuations are 
universal since our simple argument of section [5] showed that all quantum 
effects have to be of order /h. 



7.2 conductance fluctuations and speckle correlations in optics 

Here it is quite useful to compare the physics of electronic transport with 
similar physics in optics where one measures the fluctuations of a transmission 
coefficient. In optics the light scattered by a diffusing medium forms a speckle 
pattern on a screen, and we want to describe the speckle fluctuations (This is 
exactly a generalization of Young experiments. Two slits produce well defined 
fringes. Here the diffusing medium, like e.g. colloidal suspension, produces a 
complicated pattern called a speckle). A laser beam is sent on a diffusing 
medium along an incident direction a and the diffused intensity is measured 
along a direction b. The speckle pattern displayed on figure [TT] represents the 
intensity measured along a direction h for a fixed incident direction a. So a 
given intensity on the screen represents the transmission coefficient Tab from 
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Fig. 11. j4 typical speckle pattern. The white and noisy curve represents the angular 
dependence of the light intensity along the cut represented by the dashed line. The 
relative fluctuations are of order unity (courtesy of G. Maret). 

an incident direction a to an emergent direction b. We notice immediately that 
there are black spots, meaning that the relative fluctuations of this coefficient 
are of order 1. This is the Rayleigh law : 



In electronics however, relative fluctuations of the conductance are very small, 
of order 1/ ■ Here we want to understand why fluctuations are large in optics 
and small in electronics, namely to compare conductance fluctuations and 
fluctuations of transmission coefficient. 

A very convenient way to link the two fields of optics and electronics is to 
use the Landauer formalism, which explicitly expresses the conductance as a 
transmission coefficient. This is formalized by Landauer formula g — X)a& "^"-b 
The main difference between optics and electronics is that in optics, it is 
possible to measure each transmission coefficient Tab while in electronics the 
conductance is related by eq. psp to a sum over all incoming and outgoing 
channels. 

The average transmission is a probability. It is a sum of contributions of 
paired trajectories. Assuming that the angular (or channel) dependence of 
Tab is negligible, all the channels contribute equally to the conductance, so 
that, from ^ : 



6T^, 



ab 




(53) 



G^^M^Tab 

n 



(54) 



and the average transmission coefficient is thus equal to 
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Now we want to calculate the correlation between two transmission co- 
efficients, that is the function TabTa'b' ■ This quantity is the product of two 
average transmission coefficients plus a correlation term : 



TabTa'b' = Tab Ta'b' +STabSTa'b' ■ (56) 

The correlation term is constructed by pairing of trajectories corresponding 
to different transmission coefficients (Figure [12] ) . 





yb 



Fig. 12. Schematic representation of the two contributions to the product TabTa'b'- 
The first (a) corresponds to the product Tb T!i'b' ■ The second (b) gives a contribu- 
tion to the correlation function that we shall denote by STab STa'y ^ ■ 



One sees that there is no dephasing between the diffusive paired trajec- 
tories, expect outside the sample since a and a' (b and b') may correspond 
to different incoming (outgoing) channels. The second term is therefore of 
the form Tab f{o-,0'',b,b') where / is a short range function which vanishes 
rapidly as soon as a ^ a' or & 7^ 6'. If a and h are angular directions, / is a 
rapidly decreasing function of the angles. In a wave guide geometry, where 
the modes are quantized, / = 5aa'5bb'- This contribution to the correlation 
function 5Tab STa'b' is called Ci : 

STab STa'b'^^ ^ Tab Saa'Sbb' ■ (57) 

For a = a' and b = b', we obtain the Rayleigh law ((53|) . that is 

Tib = 2 Tag . (58) 

The amplitude of the fluctuations is of the order of the average. This explains 
why there are black spots on figure [11] 

In order to calculate the conductance fluctuations, we have to sum over 
all incoming and outgoing channels : 

_ _ o2 

(5g2 ^ STabSTa'b' - ^--'^bb' = M^Tag = ^ « 1 • (^9) 

aa'bb' aa'bb' 
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The sum is small since most of the terms are negligible. So our picture explains 
the important fluctuations of Tab but not the amplitude of the conductance 
fluctuations. This means that additional contributions to the correlation func- 
tion 6Tab6Ta'b' may have been forgotten. 

A next contribution is obtained by pairing trajectories in a different way. 
One possibility is to exchange the quantum amplitudes, and to have one 
crossing as shown in figure flSl a. This contribution is smaller by a factor 1/g , 
so that its contribution to the correlation function STabSTa'b' is small. But it 
has an angular dependence different from the previous contribution. Figure 
1131 a shows that there is a phase factor either for the incoming or the outgoing 
beam, so that instead of (|57p . we have for this second contribution, usually 
called C2, 

STabSTa'b' ' = ^Tab {Saa' +Sbb') ■ (60) 

The factor 2/3 results from an integration over the position of the quantum 
crossing [I]. A sum over all incoming and outgoing channels gives for this 
contribution : 




(Kf) (k;) 




(Kt) (K'_,) 

(a) (b) 

Fig. 13. Schematic representation of the C2 and C3 contributions to the correlation 
function STabSTa'b' ■ 1^1 and K3 are Diffuson contributions. K2 and K3 are Cooperon 
contributions. 



Sg' 



aa'hb' 



2 



Tab {Saa' + Sbb') 



i3_ 
3M 



< 1 



(61) 
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Still this contribution cannot explain the amplitude of the observed conduc- 
tance fluctuations, since it vanishes in the large M limit. Let us consider the 
next contribution shown on figure [T51 b. with two quantum crossings. We see 
that this so-called C3 contribution is smaller than the first one by a factor 



The 2/15 factor comes from the integration over the positions of the two 
quantum crossings |Tj (see next section). This term has no angular depen- 
dence, so that the summation over channels contains now terms and the 
contribution to the conductance fluctuations is 



The fluctuations are universal, independent on the strength of disorder. In 
summary, they are universal since the corresponding correlation function 
is constructed with two conductances and two quantum crossings, giving 
X 1/g^ ~ 1. The contributions with no crossing or one crossing cancel 
because of angular dependences. The next terms with n crossings are negli- 
gible, of order l/g"~^. In optics when one considers a speckle pattern, the 
first contribution is the most important, and the one crossing and two cross- 
ings contributions are very difficult to observe (They can be observed since, 
although very small, they are a different angular dependence, and also dif- 
ferent temporal or frequency dependences [1]). In electronics, only the third 
contribution with two crossings is important after summation over incoming 
and outgoing channels. In summary, X^afe has much smaller fluctuations 
than Tab- We have now a simple recipe to evaluate average quantities or 
correlations functions : each quantum crossing gives a factor l/g. 

7.3 Amplitude of the conductance fluctuations 

In order to calculate quantitatively the conductance fluctuations, their depen- 
dence on geometry or external parameters, we must analyze more precisely 
the structure of the paired trajectories in figure [T51 b. In addition to the cross- 
ings, there is a loop. And we have to integrate on the distribution of loops, 
like for the weak localisation correction. In contrast with the weak localiza- 
tion correction, this loop is formed by two crossings instead of one. So here 
for a given position of one crossing, we have to integrate on the position of 
the second crossing. Since this second crossing is necessarily along the loop 
formed by the two crossings, the integration over the position of the second 
crossing gives a volume element proportional to the length of the loop. For 
a trajectory of length vpt, the volume is vptX'^^. Moreover, P{t) contains 
two-Diffusons and two-Cooperons contributions. A careful examination of the 
possible crossing and trajectories shows the only possible diagrams shown on 



1/9' ■■ 



STabSTa'b' 




(62) 




(63) 
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figure [T51 Taking into account their degeneracy [1], we obtain an expression 
which is as simple as the weak localization correction : 

5G^ = 6s^i^^j t[P,i{t) + P,,,t{t)]e-"^^^ (64) 

There is an equal contribution of loops with Diffusons or Cooperons. In a 
magnetic field, the Cooperon contribution is suppressed so that the variance 
is reduced by a factor 2, as seen on fig. [TOlc. This happens for the same 
magnetic field as the destruction of the weak localisation correction, that 
is for a flux quantum through the system or through L^. 

Incoherent processes not only destroy the Cooperon contribution but also 
the Diffuson contribution. This could appear surprising since we have seen 
that this Diffuson contribution corresponds to classical diffusion, and there- 
fore should not be phase sensitive. However, this Diffuson contribution is not 
really the classical contribution, since it is constructed by pairing trajecto- 
ries corresponding to different realizations of the system. If there is a phase 
breaking event, it affects equally one amplitude and its complex conjugate. 
But here, the phase breaking event may affect differently the two amplitudes 
since they correspond to different systems. 

Like we have done above for the weak localization correction (|34|) , we can 
now evaluate quite easily SG^ given by ([M)) for different geometries from the 
corresponding expression of the return probability P{t). Let us do it for a 
quasi-lc? mesoscopic wire. P{i) is given by 

P(q,t)=5]e-^^^' 

where q — mi/L is quantized by the Dirichlet boundary conditions corre- 
sponding to a perfectly connected wire (no zero mode, see section [6.4^ . In- 
serting P{t) in ([M]) . we get 

V=6^^ = A . (65) 

Tl>0 

Let us remark that the choice of the boundary conditions is very important. 
If the system were closed or poorly connected, corresponding to Neumann 
boundary conditions, then the contribution of the zero mode would lead to 
a divergence, or at least non universality of the ffuctuations. 

In the macroscopic limit, when 3> i, we can treat the system as 
infinite, replace P{t) by its dependence (|20|) for an infinite system {td /47rt)'*/^ 
and multiply by the exponential decay e"*/"^* in the integral. Then we obtain : 
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We are not surprised by this result. It is exactly the one anticipated from our 
simple argument (|5ip treating the fluctuations as due to incoherent contri- 
butions of correlated regions of size L^. 

Finally, it is quite easy to notice, that correlations functions of n conduc- 
tances imply 2n— 2 crossings. Therefore the n*'' cumulant of the conductance 
distribution is of order 5"/(7^"~^ c>c l/g"'~'^. It vanishes for n > 2 in the metal- 
lic limit g ^ oo, so that the conductance distribution is indeed gaussian. 



8 Diffusion on graphs and spectral determinant 

The calculation of the weak localization correction or of the conductance fluc- 
tuations, as well as of other thermodynamical quantities like orbital magnetic 
susceptibility [TJ HH] can be extended to the case of any structure - called 
a network - made of quasi- one-dimensional diffusive wires. First, we note 
that the quantities of interest, like the weak localization correction ([34]) or 
the conductance fluctuations (l64l). have the same structure : 



t°'P{t)e-^*dt , (67) 

where 7 = l/r^. From ((38|) . the time integral of P{t) can be straightforwardly 
written in terms of a quantity called the spectral determinant 5(7) : 

/ diP(i)e-^* = ^-— = — ln^(7) , (68) 
Jo „ E„ +7 d-f 

where S'(7) is, within a multiplicative constant independent of 7 : 

Sh)=l[(^ + En) , (69) 

n 

En being the eigenvalues of the diffusion equation p7p . Using standard prop- 
erties of Laplace transforms, the above time integrals can be rewritten in 
terms of the spectral determinant, so that the weak localization to the con- 
ductivity and the conductivity fluctuations respectively read : 

^^ = -2s^^ ^\nS,nth) (70) 
h il 07 



(<5a2) = - 6^2^^ 1^ [\uS,i{^) + \nS,nt{.l)] , (71) 

where Sd and Smt are the spectral determinants associated respectively to 
the diffusion equation for the Diffuson and the Cooperon. f2 is the total vol- 
ume of the system. These expressions are quite general, strictly equivalent to 
expressions p4|64p . Their interest is that, on a network, the spectral determi- 
nant takes a very simple form. By solving the diffusion equation (|35p on each 
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link, and then imposing KirchofF type conditions on the nodes of the graph 
with TV nodes, the problem can be reduced to the solution of a system of N 
linear equations relating the eigenvalues at the N nodes. Let us introduce the 
N X N matrix M : 

M„„ = V coth(r/„/3) , M^p = - (72) 
^ smii rja/3 

The sum J^p extends to all the nodes /3 connected to the node a; la/3 is 
the length of the link between a and p. rjajS = lap/L^. The off-diagonal 
coefficient Map is non zero only if there is a link connecting the nodes a and 

9af3 = (47r/(/)o) Ja ^-dl is the circulation of the vector potential between 
a and f3. Nb is the number of links in the graph. It can then be shown that 
the spectral determinant takes the very convenient form |TJ 118) : 

= ( 1" ) n sinh?/a/3 det M (73) 

^ (a/3) 

Lq is an arbitrary length independent of 7 (or Lfj,). We have thus transformed 
the spectral determinant which is an infinite product in a finite product 
related to detM. Using ([75)l . mesoscopic quantities (|70l7ip can be easily 
predicted for any geometry of diffusive networks [TJ [TB] . 



9 Interaction effects 

Until now we have not considered the role of electron-electron interactions. 
They turn out to give small corrections to transport quantities like the av- 
erage conductivity, but they play an important role to understand thermo- 
dynamic properties like persistent currents [T]. Moreover, until now we have 
introduced by hand a phase coherence time (or length i^). We wish now 
to understand the microscopic origin for the loss of quantum coherence. This 
phase coherence is limited by the interactions with other degrees of freedom, 
in particular other electrons through their mutual interaction. We want to 
discuss now how e-e interactions break phase coherence. 

On one hand, interaction effects can be considered as negligible. We know 
from Landau theory of Fermi liquids that in an interacting electron gas, 
free particles have simply to be replaced by quasiparticles which are dressed 
objects, screened by the could of other electrons. These quasiparticles have a 
long lifetime which diverges when approaching the Fermi level. From Landau, 
we know that 
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where e = E — Ep is the energy of the quasiparticle measured from the Fermi 
leveL On the other hand, in a disordered metal, electrons move diffusively, 
that is very slowly. They spend long time close to each other. Qualitatively, we 
can expect that diffusion somehow enhance the effect of interactions. We may 
ask if and how expression (174]) is changed because of the diffusion. Moreover, 
since the diffusive motion is dimensionality dependent, the modified lifetime 
should also depend on this dimensionality. 

Interaction between electrons is expected to have two major effets : 

- First, each electron is not only sensitive to the disordered potential but 
also to the fluctuations of the electronic density due to other electrons. 
This additional fluctuating potential modifies the position of the energy 
levels, especially near the Fermi level. So we may expect a modification in 
the repartition of the energy levels, that is a change of the density of states 
near the Fermi level. We shall show that the density of states exhibits a 
decrease, the so-called " Altshuler-Aronov" anomaly. This reduction of the 
density of states is accompanied by a reduction of the conductivity. 

- Since the e-e interaction is a inelastic process, each quasiparticle has a fi- 
nite lifetime which limits the phase coherent properties like weak localiza- 
tion, since the coherence between time-reversed trajectories in necessarily 
limited by this lifetime. 



9.1 Screening 

In order to describe interaction effects, let us start with a few reminders 
about the screening of interaction. The bare Coulomb interaction potential 
is Uq{R) = e^/i?, that is in 3d, Uo{q) — 47re^/g^. The screened Coulomb 
interaction U (q) is given by 

^ + Xo{q)Ua{q) 

where, in the small q limit (we are interested in the diffusive regime where 
qle 1), the susceptibility Xoio) is the Pauli susceptibility, that is the density 
of states 2po. Therefore, the screened interaction is given by 

U{l) = ^—^ , U{R)^-e--^ , (76) 
q ] K IX 

where the Thomas-Fermi vector k (inverse screening length) is 

= Sttb^po ■ (77) 



In the diffusive limit qle 1, the screened interaction can by approximated 

by 
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U{q) = ^ = - = ^ , that is UiR) = . (78) 

Xo 2po 2/90 

On the scale of diffusion, the screened interaction can be considered as a local 
interaction. 

However, it turns out that screening is not instantaneous since electrons 
have to diffuse to screen a local charge. Therefore the interaction is actually 
frequency dependent. This dynamical screening is described by the frequency 
dependent susceptibility Xoi^j'^) which accounts for the dynamical charge 
reorganization : 

Xo(g,^) = 2po ■ ■ (79) 

—lUJ + Uq^ 

Therefore eq. (|78|) becomes 



U{q,uj) = - = - . 80 

Xo 2/90 Dq^ 

9.2 Density of states anomaly 

As we have done above, we shall avoid technicalities of diagrammatic theory, 
and try to get the important results from qualitative arguments. Although 
we have not elaborated on the theory of Green's functions, let me remind 
you at least that the Green's function is related to the density of states by 

pie) = ^ JlmG,iro,ro)dr„ . (81) 

As we have seen in section [31 G(ro,To) is the sum of contributions from 
all closed trajectories from ro to Tq (figure [TH a). All these amplitudes have 
different different and random phases and their contribution cancels in aver- 
age. What remains is the contribution of short trajectories, giving an average 
density of states p{e) — 2po- 

How can we effects of diffusion + interactions appear on the density of 
states? The non-interacting density of states (or Green's function) is a single 
electron property, and therefore involve single trajectories. In the presence 
of electron-electron interaction, each electron trajectory can be paired with 
the trajectory of a second electron, with which it interacts. Then by pairing 
these two trajectories, we can construct a Diffuson. More precisely, we pair 
an amplitude corresponding to one electron to the conjugate amplitude corre- 
sponding to another electron. Their interaction is represented by a wiggiy line 
in Figure [T4lb. There are actually two possible contributions, depending on 
the position of the interaction line. They are nothing but the Hartree and the 
exchange (Fock) contributions. We can conveniently separate the " diagrams" 
in three different parts : 

- A short range part close to the point Tq. It ressembles somehow to a 
quantum crossing, with a dephasing between the three trajectories. 
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- two long ranged Diffusons, 

- an interaction region between r and r' . 

So, we can easily construct the structure of this additional contribution to 
the density of states : 




Fig. 14. a) Diagram for the non-interacting density of states. After disorder av- 
eraging, the contribution of diffusive trajectories vanishes because of their random 
phase. Hartree (b) and exchange (c Fock) diagrams for the density of states anomaly. 
These contributions survive disorder averaging since they contain paired trajecto- 
ries. Upper : schematic representation exhibiting clearly the three regions, a short 
range region, diffusive trajectories and the interaction region. Bottom : usual dia- 
grammatic representation. These two representations are equivalent. 




Xp is a dimensionless parameter which describes the strength of the inter- 
action. It contains both the Hartree and exchange contributions. The 1/g 
reduction comes from the quantum crossing. And the distribution of loops 
formed by the paired trajectories is given by the return probability P{t). For 
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an energy e, the two amplitudes must be taken at different energies lo and 
oj — e. Tlierefore the two trajectories can stay in phase only during a time 
A ~ h/ 1, so that the upper cutoff in the integral is ft/e. 

Here our aim is simply to present the structure of the result without 
entering into details. Another qualitative derivation is given in [19], a detailed 
discussion is proposed in [1] , and the original calculation is done in [20] . We 
see that the amplitude and the structure of this correction to the density 
of states looks very similar to that of the weak localization correction ([34| , 
except that the upper cutoff is not but fi/e. Unlike for the weak localization, 
this correction depends on the classical return probability, so that it is not 
suppressed by a magnetic field [21! . 

The form ([82]) is approximate. A more sophisticated calculation replaces 
the upper cutoff by a Fourier transform. Moreover since g oc pqQ/td, eq. 
becomes 

A f^" 

^P(e) = 77/ P (t) cos et dt . (83) 

From the expression ((20| of P{t), we obtain the energy dependence of the 
density of states anomaly : 



Spie) oc { 



' quasi — Id 



In-r (84) 



1 1 



where the characteristic length is — ^jKDjt. The structure of this 
correction related to P(t) is similar to the weak localization correction ((39| . 
except that the characteristic length has been replaced by L^. This re- 
duction of the density of states can be observed experimentally by tunnel 
measurements, since it reflects as an anomaly in the voltage dependence of 
the tunnel conductance Gt- At zero temperature, the relative correction to 
the tunnel conductance is given by 

Gt Pa 

At finite temperature, it is not difficult to generalize the expression (|83p 
of the density of states anomaly as 

5p{e,T) = - / f{e-u;)5p{u;)duj , (86) 



where /' is the derivative of the Fermi function. After a Fourier transform, 
we find 
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Me,r)--^y RT{t)P{t) cos et dt , (87) 



where the thermal function Rrit) is given by Rrit) — nTt/ smhTrTt. The 
temperature dependence of the tunnel conductance anomaly, also called zero- 
bias anomaly is 



SGt{V,T) _ 1 
Gt Po 

or, after a Fourier transform : 

6Gt{V,T) Xp 



Gt 2TTpon Jq 



dp{e,T)f{e-eV)de , (88) 



R^{t)P{t) cos eVtdt , (89) 



This correction has been measured for various systems with different dimen- 
sionalities and the l/VV, In V, and a/F predicted by (|84l89p respectively in 
1, 2 and 3 dimensions have been observed [22j . 

9.3 Correction to the conductivity 

Taking into account the interaction between electrons leads also to a reduc- 
tion of the conductivity. Without going into the details of the calculations, 
we can argue that this reduction is a consequence of the correction to the 
density of states. Both effects result from the scattering of an electron by 
the charge fluctuations induced by disorder. The temperature dependence 
of the conductivity cr{T) is related to its energy dependence at T = OK by 
a{T) = — J f {e)(T{e)de, where f'{e) is the derivative of the Fermi function. 
Since the conductivity is proportional to the density of states (Einstein rela- 
tion), we expect that the density of states anomaly leads to a correction of 
the conductivity given by 

S.iT)_ f^jynSpie^ , (90) 



o-o J \ de J Po 



where ctq is the Drude conductivity (1221) . ^'^^ ^ static interaction, the density 
of states correction 6p{e,T) is given by (|87p . Upon Fourier transforming, we 
have 



6aiT) - -A. j RUt)Pit)dt (91) 

where Ao- is a parameter which depends on the interaction |23| . Using the ex- 
pression ((20)) of P{t), we obtain the temperature dependence of the correction 
to the conductance 
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T 



quasi 



Id 



5g{T) cx -A,, < In 



(92) 



where the thermal length Lt is defined by Lt = yJhD/T. In 2d, the temper- 
ature dependence is logarithmic like the weak localization correction. Unlike 
the weak localization correction, this correction to the conductivity is not 
sensitive to a magnetic field. Therefore both corrections can be separated 
experimentally by the application of a magnetic field. 



9.4 Lifetime of quasiparticle 

Consider a Fermi sea and inject a quasiparticle in a state \a) with energy 
above the Fermi sea. It interacts with another particle (I7), £7) and the final 
state consists in two quasiparticles (|/3), e^) and e^) above the Fermi sea 
(figure [T5)). Energy conservation implies eo. + — + e^. The lifetime of 
the state \a) is given by the Fermi golden rule and it is related to the matrix 
element of the interaction : 



I oe> • 

s 




Fig. 15. A quasiparticle in a state \a) of energy ta = e interacts with another 
quasiparticle I7) of energy — e' m the Fermi sea. The final state is made of two 
quasiparticles above the Fermi sea and one hole. 



= 2^s^|(a7|i7|/?(5)|25(e, + e7-e^-e5) , (93) 

with the constraint that < 0, e/3 > and es > 0. If the matrix element 
does not depend on energies, we see immediately that the Landau depen- 
dence comes simply from energy constraints : basically three final states can 
be chosen in a range of energy e, with the constraint of energy conservation, 
whence the dependence. Indeed the matrix element of the interaction is 
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energy independent in the ballistic case, but this is to true anymore in the 
diffusive case. The goal of the following pages is to calculate the energy de- 
pendence of the litetime in the diffusive regime. Since we do not specify a 
given state eq. ([M)) can be rewritten for the lifetime at a given energy e 

111 Otti 

(94) 

where — pQf2 is the density of states. By introducing two energy integrals, 
it can be rewritten in the form 

diu de'Y \{aj\U\l3S)\''K^-^a)Sie'-e^)S{e-Lu-ep)S{e'+Lu^ 

(95) 

where e > 0, e' < and e — w>0, e'-|-w>0 are respectively the two 
energies of the initial states and of the final states. If the matrix element is 
energy independent, we trivially recover the dependence. Upon averaging 
over disorder, the lifetime has the form 

^ A-Kvl I ujW'^{uj)dio (96) 



T-ee(e) JO 

with 

W\u;) = ^ ^ \{a^\Um\'5{t - e„)<5(e' - e^)5{e ~ uj - e^)5(e' + u: - eg) 

(97) 

The characteristic matrix element W{uj) depends only on energy transfer uj, 
but neither on e, nor on e. 

The matrix element {a'y\U\(35) can be evaluated on the basis of eigen- 
ftinction of the non interacting particles and reads 

{a-i\U\p5) ^ [ dr^dr24>l{riWJr2)^p{ri)<t)s{r2)U^{ri-r2) , (98) 



where Uuj{r) is the dynamically screened potential. The combination of wave 
functions and S function can be rewritten in terms of Green's functions 

--lmG{r,r')^y2ra{r)Mr')S{e-e^) , (99) 

a 

so that W'^{u!) can be rewritten in the form 

W^{lo) = f dridr2dr[dr'2UUri - r2)U^{r\ - r'^) 

X ImG,(ri,r'i)ImGe-^(r'i,ri) ImGe' (ra, r'^)lrivG,,+^{r'^, r^OO) 
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where the average the the product of four Green's functions has been decou- 
pled into the product of two average values. By pairing Green's functions, 
it is possible to show that the average product contains a long-range part 
related to the probability P (see eg. [THand 1 ) 



IniGe(r, r')lmGe-uj{r', r) irpa ReP(r, r' , uj) . (101) 

We deduce 

^^(^) = 2 2n2 f dridr2dr[dr'Mri - r2)U^{r[ - r'^) 

TT VqH J 

X RePrf(ri,r;,c^)RePrf(r2,r'2,-t^) (102) 
or, upon Fourier transformation : 



q#0 



At the diffusion approximation, the dynamically screened potential is given 
by (Uni) so that 

W^{uj) = -^y ^-^—^ , (104) 



which can be expressed in terms of the return probability P(t) 

1 1 

T4^2(^)^_^_ / p(t)smujtdt . (105) 



^0 

Finally, the electronic lifetime ([96|) is given by 

1 2 f°° P{t) . 2 ei 



sin^ -dt (106) 



Tee(e) TTZ^o Jo t 2 

For a metal of volume Q, we can identify two different regimes : 

9 e ^ Ec where Ec is the Thouless energy. This corresponds to time 
scales t td. In this case an electron described as a diffusive wave packet 
is insensitive to the boundaries and behaves as in an infinite medium where, 
according to P{t) = Q / {AirDtY/'^ . We obtain for the integral (fTUS]) : 



dcd 1 I to 



^0 

so that the electronic lifetime is equal to 



— - - A ( — 



d/2 



(e » -Be) 



(108) 
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where Ci = \f2lTp'^ C2 — l/47r'^, C3 = \/2/67r^. Such a behavior has been 
indeed observed in silver wires (d — 1) for which W'^iuS) cx oj"'^/^ and 
1/Tee(e) oc 6"'^/^, although the measured prefactor came out to be larger than 
the value predicted here pS] , 

• The limit e <C i?c, that is i S> t/j, corresponds to the ergodic regime in 
which the diffusive electronic wave packet explores all the accessible volume 
Q. Thus we would expect P(t) to be driven only by the zero mode. This is 
not so, because in expression (jl04p this mode has been removed in order to 
ensure electronic neutrality. The excitation energy e is smaller than and 
it is not possible to replace the sum (|104p by an integral. In this limit, we 
obtain 

-^VO-^l^f . (109) 

where the coefficient is defined by the series 

E („2+i+„2)2 ■ (110) 

The ratio E^j A is the dimensionless conductance g. For uj <C E^-, the charac- 
teristic matrix element of the interaction is thus energy independent and of 
order Ajg. The inverse lifetime in this case is [53] 



^ ( ^ ) (e « E,) . (Ill) 



= (e) 2^5 \E, 



9.5 Quasiparticle lifetime at finite temperature 

In Landau theory it is well-known that the quasiparticle lifetime at zero 
energy e = and finite temperature T is simply obtained by replacing e by 
T, so that it varies as T^. The diffusive case is more subtle. It turns out 
that in this case we cannot simply substitute e by T. This is wrong in low 
dimension. Let us see why. 

What is changed at finite temperature? The Fermi golden rule is modified 
to account for Fermi factors. The condition of filled or empty states has to 
be replaced by Fermi factors and eq. (|95p generalizes as 

Ani^^ / duj de'F{e,e',uj)W^{Lj) (112) 



Tee(e,T) 

where F{e,e',ijj) is a combination of Fermi factors fe — l/{e^'^ + 1) : 

F{e,e',L0) = /,,(! - /e-^)(l - + (1 - fe')fe-ufe'+u ■ (US) 
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The first term in this expression is larger when e > 0. It describes the decay 
of an electron-Hke state above the Fermi level. The second term dominates 
when e < and describes the decay of a hole-like state into the Fermi sea. 
For e = 0, both terms are equal. Integrating upon e', we obtain 

1 f°° e'^'^ + 1 

:^^^4..^,lj..wH.)f.^.-^^ . (114) 

This lifetime can also be obtained from the imaginary part of the self-energy 
of a quasiparticle in the presence of a screened interaction |27| . At zero tem- 
perature, we recover the result ()108p . 



9.6 Quasiparticle lifetime at the Fermi level 

We now consider more specifically the lifetime of a quasiparticle at the Fermi 
level (e — 0) and at finite temperature. Physical properties such as conduc- 
tance are expressed in terms of single particle states at the Fermi level. It is 
thus essential to understand the range of validity of the description in terms 
of independent quasiparticles. From relation (|114p . we have [28] 

Stti'^ / dujW^iuj) . . (115) 



Tee(r) "Jo ' 'Sinh/3CJ 

For the diffusion in free space, the matrix element W'^{uj) is proportional to 
w'*/^ (relation [Ml), so that 



1 ndCfi f°° dio ^ '■' ^ ''/^ 



Tee(T) 2vq Jq ujsmh/3uj \Ec 

Therefore, in three dimensions, we have 



(116) 



1 V2 r°° duj fojV'^ T F du I UJ 



, (117) 



Tee(T) in'^i/o Jo UJ sinh Puj \Ec J i^o Jo uJ^ \E, 

that is 

-j^ / \ 3/2 



up to a numerical factor. Note that the exponent of the power law is the same 
as the exponent for the energy dependence of the lifetime at zero temperature 
P08p . This result follows at once if we notice that relevant processes in the 
quasiparticle relaxation described by ujW'^{uj) are those for which the energy 
transfer uj is of order T. 

It would be tempting to generalize this result to any dimension and to 
conclude that l/Tee{T) oc T'^^'^. This is not correct for d < 2. In this case, the 
contribution of e-e processes with low energy transfer ~ dominates and 
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leads to a divergence in the integral (jll6p . In order to cure this divergence, it 
is worth noticing that T^eiT) represents precisely the lifetime of an eigenstate, 
so that the energy transfer ld cannot be defined with an accuracy better than 
1/Tee- Consequently, there is no energy transfer smaller than l/ree(T), so 
that the integral l|116p needs to be cut off self-consistently for u) smaller than 
^/Tee{T). For d< 2, we thus obtain a self-consistent relation for Tee ■ 

d/2 rp .T J, , / , , \ d/2 



UJ 



Tee(T) Vq J-^/^^^ LOSl-ahfiuj \Ee) J l/r,, Oj'^ \Ee ) ^^^^^ 

where the thermal factor has been replaced by a cutoff at ^ T. In two 
dimensions, l/Tee(T) is proportional to the temperature (within logarithmic 
corrections) : 

-J—^A^ln^ (d = 2) . (120) 

Tee[T) Ee A 

In one dimension, and since Tree S> 1, the integral becomes proportional to 
i/f^ so that the self-consistent relation leads to 

1 f Er\^^^ f T 



9.7 Phase coherence 

The time ()115p has been defined as the lifetime of a quasiparticle, generalizing 
the notion introduced by Landau to the case of a diffusive system in d dimen- 
sions. We should now evaluate the phase coherence time T0(T) which limits 
coherent effects like the weak localization correction ((34|) . This time can be 
interpreted as the lifetime of the Cooperon. Its derivation consists in calculat- 
ing directly the dephasing (e**^*)) resulting from electron-electron interaction 
and accumulated between time reversed conjugated multiple scattering se- 
quences. To that purpose, the interaction between electrons is replaced by an 
effective interaction which describes the coupling of a single electron to the 
electromagnetic field created by the other electrons. This calculation [l2l[29] 
is not developed here; see [1] for a detailed derivation. 

An alternative and qualitative approach is to consider that phase coher- 
ence is limited by the lifetime of quasiparticles. Since the multiple scattering 
trajectories that are paired in the Cooperon are defined for a given energy 
state, they cannot interfere for times larger than Tee{T). This results in an 
irreversible dephasing between the trajectories and thus a loss of phase co- 
herence. It is therefore natural to assume that 



(122) 
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Indeed, the temperature dependences predicted in equations (|118ll20ll2ip 
have been confirmed experimentally, in all dimensions by weak localization 
measurements. 

It turns out that not only these two characteristic times Tee{T) and t^{T) 
are equal (within a numerical factor), but also that the two processes, quasi- 
particles relaxation and phase relaxation, are very similar. Finally let us re- 
mark that the introduction of the low-energy cutoff in (|119p may appear as a 
handwaving and artificial way to handle the low energy divergence. We have 
shown recently that the profound reason for this divergence is that, for d < 2, 
relaxation of quasiparticles as well as the phase relaxation are not exponential 

m- 

Acknowledgments - Many points of view presented in these lectures 
have been developed during a long collaboration with Eric Akkermans, and 
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